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Variable resolution methods have become frontline CFD tools, but in order to take 
full advantage of this promising new technology, more formal theoretical development is 
desirable. Two general classes of variable resolution methods can be identified: hybrid or 
zonal methods in which RANS and LES models are solved in different flow regions, and 
bridging or seamless models which interpolate smoothly between RANS and LES. This 
paper considers the formulation of bridging methods using methods of two-point closure 
theory. The fundamental problem is to derive a subgrid two-equation model. We compare 
and reconcile two different approaches to this goal: the Partially Integrated Transport 
Model, and the Partially Averaged Navier-Stokes method. 


N omenclature 


k turbulent kinetic energy 

k r resolved turbulent kinetic energy 

k c unresolved turbulent kinetic energy 

t integral length scale 

E energy spectrum 

T energy transfer rate 

T energy flux 

F c modified energy flux defined in Eq. (29) 

P total energy production rate 

P r energy production rate due to resolved scales 

P c energy production rate due to unresolved scales 

Sij mean strain rate tensor 

Ui mean velocity 

S 2 second invariant of strain rate tensor 

C a constant 

e dissipation rate 

e r resolved dissipation rate 

e c unresolved dissipation rate 

v kinematic viscosity 

II energy production spectrum 

9 time scale of turbulence 

rj turbulent frequency 

k wavenumber 

Ki partition wavenumber 

k c filter wavenumber in LES 

Kd Kolmogorov scale 
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I. Introduction 


In the context of turbulence modeling, a variable resolution method is one that resolves more or less of 
the turbulent fluctuations in different flow regions. This can be accomplished most simply by a zonal or 
hybrid method, in which RANS equations, which model all fluctuations, are solved in some flow regions, 
and LES, which models only small-scale fluctuations, are solved elsewhere. An alternative is a bridging (or 
seamless , or continuous) model, which interpolates smoothly between RANS and LES. Variable resolution 
methods for turbulent flows attempt to fill the gap between industrial RANS methods, which may sometimes 
be too coarsely resolved to be useful, and LES, which might resolve unnecessary details, thereby incurring 
excessive computational cost. At present, however, variable resolution methods cannot be said to have 
attained the level of theoretical rigor of either RANS or LES. To contribute to a more sound understanding 
of these methods, this paper outlines the introduction of spectral closure theories as a foundation for variable 
resolution methods. 

It may be helpful to begin by considering what an ideal RANS/LES hybrid model would be. Suppose that 
we knew that some hypothetical RANS model was a limit of LES under well-defined analytical conditions. 
This would entail the possibility of statistical comparisons at the level at least of two-point correlations, so 
that we could state with confidence when RANS can replace LES; we would also know a priori conditions 
under which the RANS solution will become unrealistic and must be replaced by LES. Suppose further that 
the transition between RANS and LES can be gradual enough that overlap regions exist in which both 
models are valid, so that a smooth transition is possible, rather than a transition through a fixed sharp 
interface. Then it would be possible, at least in principle, to solve the RANS model wherever it is valid and 
the LES model elsewhere. This would be a zonal method in which LES and RANS are smoothly connected 
though an overlap region, the location of which is determined as part of the numerical solution. 

Since these requirements might seem somewhat unrealistic, we describe a deterministic analog, in which 
a RANS model is connected, not to LES or DNS, but to a spectral closure. 1 Like the hybrid described 
above, this hybrid connects a statistically complex, more realistic description (the spectral closure) to a 
simpler abridged one (the RANS model). In this case, a multiple-scale perturbation analysis demonstrates 
that when 

e<e 2 /fc and £ < k 1/2 (1) 

where e is the dissipation rate, k is the turbulent kinetic energy, and £ is a length characteristic of the large 
scales, there is an approximate solution of the spectral closure of the form E(n,t) = E(n;e(t),£(t)) where 
E(n,t ) is the time-dependent energy spectrum, and where e and £ satisfy a RANS-type two-equation model 
obtained as part of the perturbation scheme. This model is, however, not applicable if the inequalities in 
Eq. (1) are not satisfied; in this case, the complete spectral closure must be solved. Since the conditions 
for an ideal hybrid are satisfied, a (temporal) hybrid model could be constructed that combines the spectral 
closure and the two-equation model, so that the complex spectral closure is solved only when required by 
the dynamics. 

This description of a hypothetical hybrid is intended only to motivate the theoretical possibility of such a 
model. Of course, in practice, or at least in current practice, the requirements of the ideal hybrid model are 
far from being satisfied: we do not know precise conditions under which LES can be replaced by RANS, nor 
are there known a priori conditions under which the RANS description must become inadequate. Instead, 
the adequacy of RANS is assessed after the fact, by comparison with numerical or experimental data. There 
is no connection between RANS and even partial statistical information like two-point correlations; thus, 
the more complete agreement at a refined statistical level needed for a meaningful hybrid method is not 
available. Interface conditions are typically prescribed according to various empirical heuristics: there is 
neither theorical guidance, nor even any general concensus, on how to switch between RANS and LES, or 
on how to treat the buffer region between them. All of these limitations mean that the validity of hybrid 
models cannot be insured in advance and problems must be fixed as they arise. 

Bridging models offer a practical way to overcome at least some of these limitations. In the simplest terms, 
the idea is that RANS can be connected to LES or DNS by reducing the RANS eddy viscosity according to 
grid size or some other numerical consideration. 2 The expectation is that reducing the eddy viscosity can 
“free up” fluctuations since in the limit of zero eddy viscosity and, of course, sufficient numerical resolution, 
DNS is certainly recovered. 

Although this idea is attractive, later work 3 showed that better control over the consistency of the RANS 
and LES models is necessary: for example, it should be possible to demonstrate that the total energy is 
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independent of the resolution. Such requirements impose much tighter connection between the RANS and 
subgrid model formulation than is possible by a simple damping of the eddy viscosity. Yet this requirement 
poses a difficulty, because the most popular subgrid formulation, Smagorinsky-based LES, is in RANS 
terminology a zero-equation model; thus, it seems difficult to connect such subgrid models smoothly, say to 
a two-equation RANS model. 

The simplest way to enforce consistency between RANS and LES is to develop a subgrid two-equation 
model that reduces to a standard RANS model in the limit that all fluctuations are modeled. Two current 
systematic approaches to this goal are the Partially Integrated Transport Model (PITM) of Schiestel 4 and 
the Partially Averaged Navier-Stokes (PANS) of Girimaji. 3 The purpose of this paper is to clarify the 
formulation of PITM with particular reference to the role of spectral closure, which appears to be a natural 
theoretical framework for this type of modeling. The arguments and conclusions of PANS will be compared 
with PITM as appropriate. Our immediate concern is to reconcile the derivations of the subgrid dissipation 
rate equation, because the same conclusion is reached in both PITM and PANS, yet the arguments are 
quite different. We conclude with some observations on the role of spectral closure and of self-similarity in 
formulating bridging models. 


II. Background from Spectral Closure Theory 


PITM originates in Schiestel’s multiple-scale model, 5 which we briefly summarize. The foundation is 
spectral closure , the Fourier form of two-point closure for homogeneous turbulence. In the usual formula- 
tion, the basic unknown is the energy spectrum E(n,t), the Fourier transform of the two-point correlation 
function, under the kinematic assumption of isotropy of fluctuations. Important single-point quantities can 
be expressed as moments of the energy spectrum, for example, the kinetic energy and dissipation rate are 


COO coo 

k= E(k)cIk e = 2v n 2 E{n)dK 
Jo Jo 


The wavenumber variable k corresponds to the intuitive idea of a “scale of motion.” 
The spectral evolution equation 

E(k, t) = II(«;, t) — T(k, t ) — 2 uk 2 E(k, t ) 


( 2 ) 


( 3 ) 


can be derived from the Navier-Stokes equations. In Eq. (3), II is a spectrum of energy input, T is the 
nonlinear transfer, a third-order moment of the velocity field, and the last term is the dissipation spectrum. 6 
Although the time variable is appears in Eq. (3), subsequently it will be assumed rather than written 
explicitly. The spectral formulation is convenient because it can model the interaction between different 
scales of motion more or less faithfully, depending on what is used for T . This makes possible a systematic 
formulation of the interactions between resolved and unresolved scales of motion. 

In studies of homogeneous isotropic turbulence, it is usual to assume that the energy input is concentrated 
at some large scale, so that II(re) is nonzero only inside some interval \k — «o| < Ak. The total production 
(energy input) rate is 

coo 

P= dn ![(«) (4) 

Jo 

A spectrum for production by mean shear can be constructed by using Lumley’s 7 inertial range spectrum 
for shear stress. The result is 

n(«) = CS 2 e 1/3 K“ 7/3 (5) 

where S 2 = SijSij and Sij = dUi/dxj +dUj/dxi is the mean strain rate, XJi is the mean velocity, and C is a 
model constant. Eq. (5) is an inertial range result only; some modification near k = 0 is assumed to insure 
finiteness. The crucial difference between the large-scale forcing used in studies of homogeneous isotropic 
turbulence and production by mean shear is that in the latter case, energy production occurs at all scales, 
not only at large scales. A simple way to motivate Eq. (5) is to note that the total production due to scales 
with wavenumbers larger than k is 


cOO 

P(k) = C S 2 / dp e 1 /: V 7/3 oc S 
J K 


2 e 1 / 3 K -4 / 3 


( 6 ) 


where e 1 / 3 ^ 4 / 3 is the scale-dependent eddy viscosity corresponding to Kolmogorov inertial range scaling. 
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Energy conservation by nonlinear interaction implies 


dn T(k) = 0 (7) 

This fact permits the definition 

T(k) = |^(k) (8) 

where T is the energy flux satisfying .A(O) = .F(oo) = 0. 

It will be useful to note some elementary models for the energy flux. The simplest is the Kovaznay model 8 



T{k) = Ctj(k)E(k) where rj = y/ k 3 E(k) 


(9) 


Here and subsequently, we adopt the convention that C denotes a constant, not necessarily the same one 
each time it appears. A more realistic model is the Heisenberg model 6-8 

pK pOO 

.F(k) = C / dp p 2 E(p) / dp E(p)9(p) where 0(p) = 1/ tJp 3 E[p) (10) 

Jo J K 

The Kovaznay model describes energy transfer as a stepwise cascade in wavenumber space; the Heisenberg 
model treats energy flux as local energy production: the product of a spectral eddy viscosity (the second 
integral on the right side of Eq. (10)) and the square of a local strain (the first integral). Some refinements 
are offered by a generalized Heisenberg model, 9 

= C ( r dp p 2 E{p) f°° dp E(p)0(p) - f dp p 4 r dp j (11) 

U o Jk Jo J k V J 

Whereas the Kovaznay and Heisenberg model treat transfer as purely “forward,” that is, from large scales 
to small scales, this model allows transfer from small to large scales, often referred to as energy backscatter. 
This feature is crucial to a correct formulation of the problem of isotropic decay. 

Integrate the spectral evolution equation over all wavenumbers k using the relations (some of which have 
already been noted) 


to obtain the energy balance 


/ dn. E(k) = k 
Jo 

poo 

/ dK T(k) = 0 

Jo 


poo 

/ d,K II(k) = P 

Jo 

poo 

/ dn 2 izk 2 E(k) = e 

Jo 


( 12 ) 


k = P — e (13) 

We note that although this equation provides an exact relation, the nonlinear transfer T ', the heart of the 
turbulence problem, plays no role in it. 


III. Formulation of PITM 

Schiestel 5 refined the energy balance of Eq. (13) by introducing partition wavenumbers «].,•■■ K n , corre- 
sponding spectral “bins” or “slices” i < n < «$, and partial energies 

ki = f dn E(k) (14) 

J 

Heuristic arguments were used to obtain coupled sets of two-equation models for the k, : and fluxes T[k J). The 
result is an alternative to numerical integration of Eq. (3) with some flux model like Eq. (10). However, this 
model is not a consequence of spectral closure, but is instead a heuristic modification specifically designed 
to connect to the usual RANS models. Our suggestion is that this multiple-scale model could be formulated 
as a perturbation theory in which the energy flux is a slowly varying function of wavenumber, so that the 
spectrum is approximately of Kolmogorov form at each wavenumber, but a global Kolmogorov spectrum 
over the entire inertial range does not necessarily exist. 
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Chaouat and Schiestel 10 applied this formalism to two-equation subgrid modeling by taking only one 
subgrid wavenumber shell. In their words, “In this framework , a cutoff wavenumber k c is introduced in the 
medium range of eddies while the wave number Kd is located at the end of the inertial range of the spectrum 
after the transfer zone. ” We assume then that their Kd is the inverse Kolmogorov scale, 


Kd 


(15) 

where v is the kinematic viscosity. With this understanding, the unresolved scales are k c 
resolved scales are k < k c \ the RANS limit is k c — > 0 and the DNS limit is k c — > oo. 
Chaouat and Schiestel 10 define the resolved energy 

< k < Kd and the 

k r = 

f dhi E(k) 
Jo 

(16) 

and subgrid (or unresolved) energy 

/»oo 


k c = 

/ dK E(k) 

(17) 


so that the sum of the resolved and unresolved energy is the total energy. The partially integrated spectral 
evolution equation contains the terms 



dre E(k) = k r 

dn T(k) = f 
Jo 


dn II(k) = P r 


dK - 7 — = J-(k c ) 

8k v ' 


dK 2 vk 2 E(k) 


€ r 


so we obtain the resolved energy equation 


k r = P r — E(k c ) 


€ r 


(18) 


(19) 


For production by forcing concentrated near ko C k c , P r = P. We expect dissipation to occur at small 
scales; then e r ~ 0, so consequently Eq. (19) becomes 


k r = P — J-(k c ) 


( 20 ) 


In the Heisenberg model Eq. (10), 


E(k c ) = C 



dp p 2 E(p) 



dp E(p)6{p) 


( 21 ) 


Although we will not solve the Heisenberg model, it suggests the usual heuristic that the subgrid scales act 
on the resolved scales as an eddy viscosity. We might model this subgrid eddy viscosity by a two-equation 
ansatz 

r°o l2 

/ dp E(p)0(p) = C— (22) 

Jk,c 

In this case, we need equations for k c and e c « e. 

The obvious starting point is to integrate the spectral evolution equation over the unresolved scales. This 
generates the terms 


so we obtain 



k c — Pc T H*c) e c 


(23) 

(24) 
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Note that the sum of Eqs. (24) and (19) is the equation for total energy, Eq. (13). For production by forcing 
concentrated near kq k c , P c = 0. For production by shear, 

Pc = ^C5 2 c 1 / 3 «- 4 /3 (25) 

is small when k c is large, but since the decay is only algebraic, the unresolved production P c can be significant. 
Again, since dissipation will occur at very small scales, we expect e c s=s e, so 

kc = P c + F(k c ) - e (26) 


A key point for Schiestel is 

that k c can be a function of time. Allowing k c = K c (t), we have 


d 

dt , 

/•OO f‘00 

/ dK E(K,t) = / dn E(K,t) — k c (t)E(K,t) 

J K c (t) J K c (t ) 

(27) 

which adds a term to Eq. (26): 

kc — Pc + F(k c ) P(ftc)k , c € c 

(28) 

Define a new quantity 

F c — F(k c ) F(j^c)kc 

(29) 

so that 

k c = P c + F(k c ) - e 

( 30 ) 

Then solve Eq. (29) for k c , 

P(k c ) - F c 
~ E(k c ) 

( 31 ) 


Here, some caution seems indicated, since what is happening is a reshuffling of definitions, not the introduc- 
tion of new relations. 

Attempts to find a dissipation rate equation from “equilibrium” considerations alone introduce the Kol- 
mogorov relation 

K = C fc372 ( 32 ) 

which is just a restatement of the Kolmogorov energy spectrum E(k) = C e 2 / 3 K -5 / 3 . Chaouat and Schiestel 10 
also use the differential form of this relation 


e 3 k k 

e 2 k k 


(33) 


Here however we come to a difficulty, because if we integrate the Kolmogorov spectrum from k c to the 
Kolmogorov scale Kd, we obtain 


kc kd — C 

Therefore, after some rearrangement, 


/ 


«d Q 

2/3 -5/3 = £ 
2 


Ce 2 / 3 (/ 


-2/3 


- K, 


-2/3' 


1 


2/3 2/3 

K c K d 


-3/2 


= c 


(k c - k d ) 3/2 


(34) 


(35) 


(recall the convention on the use of the constant C), whereas Chaouat and Schiestel 10 state a dimensionally 
consistent alternative 

K d ~ k c = (36) 

k c 


That Eqs. (35) and (36) are not equivalent is shown by the limits as Kd —■ > oo of their left sides: 


2/3 


2/3 


-3/2 


( 37 ) 


6 of 9 


American Institute of Aeronautics and Astronautics 


but 


Kd ~ K, 


00 


(38) 


At this point, a basic principle, or perhaps belief, about turbulence becomes relevant: it is that limits 
Kd — ► oo, or equivalently, in view of Eq. (15), v — » 0, are always finite. The idea is that letting v — > 0 creates 
more small scales but leaves large scales unchanged. This idea can be stated as the existence of a high Re 
“fixed point;” it permits the usual modeling heuristic that Reynolds number appears only in low Reynolds 
number models. This principle would seem to favor Eq. (35) over Eq. (36); but perhaps it is uncertain that 
this principle applies to the interval k c < k < Kd- 

Meanwhile, completing the model requires two more relations. Schiestel 5 used the Kovaznay closure 
for the (static) flux T, given by Eq. (9). This might be open to the objection that it uses both a model 
(Kovaznay’s) and its solution (the Kolmogorov spectrum), but Chaouat and Schiestel 10 do not use any 
particular closure for the flux; instead, they set 


F r = 


(39) 


to eliminate F c from the equations. 

We complete the derivation using what we think are correct equations, but following Chaouat and Schi- 
estel 10 by leaving the flux model unspecified. Return to Eq. (31) and substitute Eq. (39), 


. F(k c ) F c F^Kc') € K c . , . . 

= -%r = = - e) 


where 


k c E(k c ) 


(40) 

(41) 


F(k c ) 


- 1 


is a (possibly) scale-dependent factor, although in an inertial range, a = | is a constant. 
Next, return to Eq. (33) to obtain 

3 k c k c 3 e 3 e 2 e 2 

e = 2 Y c e + V c e = 2 Y c {Pc + Fc) - 2 V c + a k c 

This equation states a subgrid dissipation rate equation 

e = C a e ^(P c +F c )-C a J f 

rZr for 


with 


C S J> = - C s J = --<j( - 1 

el 2 e2 2 V e 


(42) 


(43) 


(44) 


Because of using Eq. (36) instead of Eq. (35), a much different result for the constants is obtained in 
Chaouat and Schiestel: 10 

3 


C\l = 


but 


(JSQ ^ 

2 KdE(Kd) 


F(k c ) 


- 1 


so that when k c = k and all of the kinetic energy is modeled, we obtain the RANS model coefficient 

_ 3 k ( F(k c ) n 

62 2 KdE(Kd) 


(45) 

(46) 

(47) 


In either case, the value = |, obtained here because of the Kolmogorov relation Eq. (33), is in reason- 
able agreement, perhaps coincidentally, with the value obtained by calibrating to homogeneous shear flow. 
However, since KdE(Kd) oc k^ 2 ^ (recall that Kd is assumed in Chaouat and Schiestel 10 to be at the end of 
the inertial range), the model constants in Eqs. (46) and (47) are divergent as Kd — < ► oo. 
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That this is unsatisfactory is of course recognized in Chaouat and Schiestel, 10 and some arguments are 
given to lead to a finite result. One possible approach is to evaluate the ratio of the last terms on the right 
sides of Eqs. (46) and (47); this (purely formal) step cancels the divergences to give 


n s 9 _ 3 , 

^e2 2 _ 

0,2 - I “ k 


(48) 


and so permits Chaouat and Schiestel 10 Eq. (46), 


c s J-^ = ^(c e2 --) 


(49) 


where now we can set 0,2 to some empirical value. We stress that although this formal manipulation 
produces a finite answer, this theory does not predict a definite value of 0,2- 

To conclude, we have obtained the subgrid dissipation rate equation in the form 


O el f (P c +F C ) -fc e i + ” 77 (0,2 - O,!)) e - 

rC c V 

0,1 y{P c + F c ~e)~ (0,2 - Cel)j 


(50) 


although the argument cannot inspire complete conhdence. It is therefore important that the same subgrid 
dissipation rate equation is obtained in PANS by a simpler, more straightforward argument: 

1. Postulate the standard equation 

e = C e \jP — 0,2 y- (51) 

k k 

k k 

2. Remark that in any self-similar flow, k c /k is constant, so — = — . Then 


P-e= -(P c + F c -e) 
kr 


(52) 


The appeal to self-similarity is justified by the fact that the dissipation rate equation is calibrated to self- 
similar flows; we will return to this point shortly. 

3. Substituting in the standard equation, 


2 2 

e = C el e -+ 0 £l f (P c +F c -e)~ C, 2 ^ 

nj (C 


(53) 


in agreement with Eq. (50) after simple regrouping of terms. This alternate derivation suggests that 
although we can question some aspects of the argument of Chaouat and Schiestel, 10 incuding Eq. (36) and 
the consequent need to cancel divergent terms (although we stress that Eq. (48) is our own interpretation 
of Chaouat and Schiestel 10 ), the final result is consistent with an alternative derivation which does not raise 
these questions. 

For the bridging model to succeed, it must indeed connect RANS to LES. In the limit k c — > 0, Eq. 
(28) reduces to Eq. (13) because n c being constant, k c = 0; the production becomes the total production, 
and the dissipation rate is unchanged. The reduction of the resolved dissipation rate equation to its RANS 
counterpart is obvious from Eq. (50) on setting F c = 0, P c = P and k c = k. 

What is perhaps less trivial is the DNS limit n c Kd- In this limit, k c — > 0; substituting Eq. (49) in Eq. 
(43) shows that in this limit, we can have a steady state in which P c + F c = e. Then the entire two-equation 
subgrid model is replaced by a local equilibrium condition and the two-equation model effectively disappears 
in the limit of high resolution. 

Recall that for production by a mean shear, the subgrid production is given by Eq. (25). For very large 
k c corresponding to highly resolved LES, approaching DNS, P c oc -C F c oc k°: direct production by 

shear is dominated by local nonlinear transfer. Thus, homogeneous shear is not an appropriate calibration 
case for large k c \ consistency with a forced steady state is more plausible. This change of calibration case is 
implied by the possibility, noted previously, of a steady state in the subgrid scales in the DNS limit. This 
consistency check gives some theoretical support to the subgrid dissipation rate equation common to PITM 
and PANS. 
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IV. Comments on the Subgrid Dissipation Rate Equation 


In a self-similar flow, quantities with the same units are constant multiples of each other, 
self-similar decay, the equation 



Thus, in 
(54) 


is not a “model;” it is rather a consequence of the assumption of self-similarity. Since different self-similar 
flows can give different values of the proportionality constant, such a relation is not at all a general principle, 
and there is therefore no reason to believe that Eq. (54) can be deduced from any theory. The standard 
dissipation rate equation Eq. (51) is simply a more general statement of self-similarity that can apply to 
two self-similar flows, usually chosen as homogeneous shear and decay: data from these flows is used to 
determine the model constants. 

As derived in Eqs. (51 ) — (53) , the subgrid dissipation rate equation follows from the usual self-similarity 
assumptions and the additional observation that k c /k is constant in a self-similar flow. Since we do not 
consider that this equation can be derived from anything more basic, the elaborate explanation in Chaouat 
and Schiestel 10 may not add anything significant to the very elementary PANS argument that is based openly 
and explicitly on self-similarity alone. 

We have stressed the role of self-similarity in formulating the standard and the subgrid dissipation rate 
equations. These equations describe certain calibration cases; their extension to other problems comes with 
no guarantees. A different viewpoint 1 is that the length-scale determining equation of the two-equation 
model explicitly must describe not self-similarity or “equilibrium” but departures from them. Woodruff 
and Rubinstein 1 propose a two-equation model based on the classical Heisenberg model Eq. (10) that is 
consistent with this idea. 


V. Conclusions 

The spectral formulation of PITM is a useful and convenient way to discuss variable-resolution modeling. 
Deriving the subgrid energy equation is straightforward; the new element in the role of nonlinear transfer in 
addition to direct energy production by shear. In PANS, the subgrid dissipation rate equation is formulated 
in terms of self-similarity (or ‘fixed-point’ analysis); the agreement with the PITM formulation, obtained by 
an apparently different argument, perhaps illustrates the extent to which the assumption of a (self-similar) 
Kolmogorov spectrum limits the range of possible end results. An interesting new feature of the subgrid 
dissipation rate equation is that one calibration case changes from homogeneous shear to steady-state forced 
turbulence as the resolution increases. Finally, the subgrid two-equation becomes nugatory in the DNS limit; 
this fact insures the consistency of the bridging model with this limit. 
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